Interlayer c-axis transport in the normal state of cuprates 
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A theoretical model of c-axis transport properties in cuprates is proposed. Inter-plane and 
in-plane charge fluctuations make hopping between planes incoherent and diffusive (the in-plane 
momentum is not conserved after tunneling). The non-Drude optical conductivity a c (ui) and the 
power-law temperature dependence of the dc conductivity are generically explained by the strong 
fluctuations excited in the process of tunneling. Several microscopic models of the charge fluctuation 
spectrum are considered. 



Despite the strongly two-dimensional layered structure 
of the high-temperature cuprate superconductors, fea- 
tures associated with the third dimension, perpendicular 
to the CuOi planes, may be an important ingredient in 
their superconductivity. In fact, it is well accepted that a 
certain degree of Josephson-type coupling between differ- 
ent planes is necessary to suppress the two-dimensional 
fluctuations, which will otherwise destroy the supercon- 
ducting long-range order. However, the systematic de- 
pendence of the critical temperature T c on the number 
of layers in the unit cell (together with the absence of ev- 
idence for strong fluctuations effects above T c at optimal 
doping, which suggests that these fluctuations are not the 
major reason for this systematic dependence) points al- 
most unambiguously to the conclusion that theories for- 
mulated for a single plane cannot be the whole story. 
Either hopping between planes JjJ], or Coulomb interac- 
tion between them [0, or both, is an important factor in 
raising the critical temperature (and, perhaps, in some 
cases also for lowering it, see Ref. M). In the light of this, 
the study of the c-axis optical and transport properties 
is more than just a minor diversion from the main issue. 

These c-axis optical and transport properties are very 
puzzling and anomalous ||. Most remarkable is the 
fact that the temperature dependence of the dc c-axis 
resistivity p c (T), in sharp contrast to the well-known 
linear T-dependence of the in-plane resistivity p a b(T), 
is non-universal, being described in most cases by a 
power law p c (T) ~ T 7 , where however the exponent 
7 can be anything in between approximately +1 and 
— 1. The optical conductivity a c (u) is roughly fre- 
quency independent from low frequencies up to mid- 
infrared frequencies (except in the case of some over- 
doped cuprates (Y \Ba2Cu3O7 and La2-xSr x CuC>4))] the 
numerical value is below the Mott-Ioffe-Regel minimum 
metallic conductivity. This behavior, which is dramati- 
cally different from the behavior of the in-plane resistiv- 
ity, has been christened "confinement" H. Thus, despite 
the dramatic differences in the raw data between the dif- 
ferent cuprate families, one can isolate at least two el- 
ements which may be legitimately called "universal" : a 
non-Drude optical conductivity of a magnitude below the 
Mott limit, and a power-law temperature dependence of 



the dc resistivity p c (T) (albeit with a material-specific 
exponent). In this paper we develop a framework for 
the explanation of these universalities, which will hope- 
fully shed light on the role of the inter-plane and in-plane 
Coulomb interaction as well as on the more obvious one 
of the inter-layer hopping. 

On the basis of the above experimental observations, 
we can make the assumption of incoherent transport (dif- 
fusive tunneling) in the c-direction. The inter-plane (or 
rather inter-unit cell) hopping time Thop can be estimated 
from the dc c-axis resistivity. Relating the diffusion con- 
stant D to the hopping time r^ op by D = c? 2 /r^ op , we can 
derive a model-independent relation between the conduc- 
tivity er c and the hopping time: a c = e?vid-^—-, where 
V2d is a two-dimensional density of states. Using this 
formula and the experimental values for the c-axis con- 
ductivity we can estimate the c-axis hopping time T>> op . 
The strong two-dimensionality of the electron motion be- 
comes obvious if we compare the hopping time Thop with 
the in-plane scattering time T a \,. As is well known, the 
in-plane scattering time T a b is of order of fi/kT. Direct 
comparison Q shows that for most materials the c-axis 
hopping time is much longer than the scattering time 
in the plane. These two times are comparable only for 
overdoped YBC12C113O7 and La2~ x Sr x CuOi suggesting 
a crossover to a different regime of c-axis transport; this is 
confirmed by the experimental observation of the Drude- 
like frequency dependence of the conductivity o~ c (u>) for 
these compounds. 

Many approaches have been suggested |l],|3| to describe 
the c-axis transport properties. Most of them stem from 
phcnomcnologically assumed in-plane Green functions. 
One remarkable example is a non-Fermi "Luttinger" liq- 
uid theory which explains naturally the "confinement" 
for c-axis motion in the normal state Q. Many other 
theories are essentially based on the Fermi liquid theory 
modified by strong correlations . 

We take a quite different approach to the problem. Al- 
though in our approach the in-plane motion (expressed 
by in-plane Green functions) is undeniably important, we 
show that most of the c-axis properties can be qualita- 
tively understood on the basis of knowledge of the spec- 
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trum of in-plane and inter-plane charge density fluctua- 
tions excited during the process of the inter-plane tun- 
neling. The spectrum of charge fluctuations can be 
directly measured experimentally by optical reflectivity 
and electron-energy- loss spectroscopy (EELS). 

The essential physical picture of our approach is that 
the c-axis tunneling is strongly suppressed by charge fluc- 
tuations excited in the process of tunneling B. This 
anomaly (the so called Coulomb blockade) is widely ob- 
served in many strongly correlated and mesoscopic sys- 
tems. Other examples of this class of phenomena are 
orthogonality catastrophes and a zero-bias anomaly in 
diffusive systems. The ubiquity of the Coulomb blockade 
phenomena (static or dynamic) in correlated systems in- 
dicates that the anomalous c-axis transport properties 
may be merely a consequence of strong correlations in 
the cuprates. The necessary condition for the appear- 
ance of the Coulomb blockade phenomenon is the strong 
effective coupling of the tunneling electron to the collec- 
tive excitations of the liquid. In fact, we can think about 
the Coulomb blockade as a "high-energy" phenomenon of 
order of Coulomb energy per electron and independent of 
the low-energy quasiparticle spectral properties, be they 
fermi- or non-fermi liquid. In other words, the tunneling 
electron couples to the excitations in the broad range of 
frequencies from low up to high frequencies. The non- 
fermi liquid property (property not present for a three- 
dimensional Fermi liquid or for a electron gas in the RPA 
approximation) which is responsible for the anomalous c- 
axis properties is simply the large density of "detuning" 
charge fluctuations (as observed by optical and Raman 
spectroscopy) over a broad frequency range up to mid- 
infrared frequencies. 

Empirically, an overview of the in-plane conductivity 
and the c-axis conductivity in various families of cuprates 
does not reveal any obvious correlation between their 
temperature dependencies. On the other hand, the per- 
turbative diagrammatic expression (assuming the equiv- 
alent Green's functions in each plane and uncorrelated 
impurities) for the in-plane conductivity and out-of-plane 
conductivity H are equivalent up to the vertex func- 
tions. Therefore the difference in the temperature and 
frequency dependence of the in- and out- plane conduc- 
tivities must come exclusively from the interplane tun- 
neling probability. For instance, the notion of the "two- 
dimensional Luttinger" liquid is not sufficient by itself to 
explain the difference between the in-plane and out-of- 
plane resistivities 0. 

One particular approach to explaining the difference 
between the resistivities along the different directions is 
based on the highly anisotropic form of the tunneling 
matrix element t±(k x ,k y ) as a function of the in-plane 
momentum. Several authors || developed a phenomeno- 
logical approach (assuming as well a strong anisotropy of 



quasiparticle lifetimes and the density of states around 
the fermi surface) which seem to fit successfully the ex- 
perimental data. Two remarks are in order. First, this 
approach assumes that the tunneling conserves the in- 
plane momentum. This may be the case in certain sit- 
uations (possibly, in the superconducting state and in 
the overdoped regime), but in general this assumption 
deserves a close scrutiny by experiment and theory. In 
fact, we will argue that the in-plane momentum is not 
conserved when r/ lop 3> r a b- Second, the anisotropy of 
the tunneling matrix element t±(k x ,k y ) is different for 
some variations of cuprates (8j, and , in general, it can 
be doping dependent. Thus in some cases, t^_(k X: ky) 
may not vanish along the diagonals of the Brillouin zone 
making the contribution of the diagonal quasiparticles 
(k x = iky) to the c-axis conductivity non- vanishing con- 
trary to the assumptions of references § . 

In the rest of the paper, we begin by generalizing the 
standard tunneling formalism, introducing a non-trivial 
tunneling probability which accounts for the inelastic and 
elastic (momentum scattering) processes. Then, we cal- 
culate this tunneling probability from assumed spectra of 
"the detuning fluctuations" . After that, we calculate the 
experimentally measured optical conductivity <t c (lu,T) 
and the tunneling conductance o~ c (V). Finally, we discuss 
the complex experimental situation and possible exten- 
sions of the proposed theory. In the appendix we discuss 
several microscopic models giving the spectrum of the 
charge fluctuations. 

Theory. The tunneling formalism. The tunneling of an 
electron from one plane to another plane can be consid- 
ered by using the time-dependent tunneling hamiltonian 
formalism. The "blockade effect" due to the excitation 
of the electromagnetic modes is accounted by the mod- 
ulation of the tunneling matrix element by the Coulomb 
interaction. Thus the part of the hamiltonian responsible 
for the transfer of electrons between planes is: 

H c = tx(rur2;t){at{ri)a 2 {r 2 ) +4fa)»i(ri)), C 1 ) 

ri,r 2 

where the quasi-classical tunneling matrix element 

tx(ri,rr,t) is equal to tx(n, r 2 )exp (-^ A(z,t)dz*j . 

Gauge invariance dictates the presence of the phase fac- 
tor <p(r\,t) = (e/c) JJ* 2 A(z,t)dz, where the integral is 
taken over a path connecting two points ri and r 2 on 
the different planes. Since the optimal tunneling trajec- 
tory is perpendicular to the planes (along the z axis), the 
tunneling matrix element can be written as tj_(ri,r 2 ) = 

t±S(ri — r 2 )exp (^— ™ J** A(z, t)dz^j . It conserves the in- 
plane momentum. Using the tunneling hamiltonian for- 
malism P,|lC| , we get the following expression for the tun- 
neling current I(t) between two planes: 
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I(t) = ~Re 

h 2 



drdr 1 



+00 



dt'e- l ^\t ± \ 2 P{r -r',t-t')S(r -r',t-t'), P(r-r',t-t') = ^'•''e-'"^), (2) 



S(r-r',t-t') = Q(t -t') < G<{r - r' ,t - t')G>{r' - r,t' - t) - G<{r ~ r' ,t - t')G>{r' - r,t' - t) > 
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where P(r — r' ,t — t') is a phase-phase correlation func- 
tion between two planes averaged over the equilibrium 
fluctuations, and V is an applied voltage. The defini- 
tions of Green's functions and essential details of the 
derivation can be found in the Ref. j[o|. The fact that 
the hopping time r^ op is much longer than the in-plane 
scattering time T a b allows us to separate the in-plane 
propagation S(r — r',t — t') and the tunneling proba- 
bility P(r — r',t — t'). It is important to remark at 
this stage that the long-wavelength fluctuation modes 
(with the wavelength much longer than the in-plane mean 
free path) can suppress the tunneling probability with- 
out effecting the in-plane motion. The scattering by the 
short-wavelength fluctuations is accounted by the spec- 
tral properties of the in-plane propagation S(r—r',t—t'). 
The understanding of the properties of the tunneling 
probability P(r — r',t — t'), describing the effect of the 
"detuning fluctuations" , is imperative for any particular 
problem of the tunneling. The importance of this cor- 
relation function was first described in Ref. pT| . The 
novel element here is the discussion of the spatial depen- 
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dence of the tunneling probability function P(r—r', t—t'). 
The spatial dependence appears to be very important for 
many questions of the c-axis transport properties. As 
mentioned above, in previous studies of c-axis transport 
in cuprates, specific properties of the tunneling probabil- 
ity (e.g. the tunneling with or without the conservation 
of the in-plane momentum k\\) were assumed. Here we 
analyze and calculate the tunneling probability from the 
fluctuation spectrum of the electromagnetic field. We 
can call the tunneling "diffusive" if the in-plane momen- 
tum is not conserved (if the momentum is conserved, it 
can be called specular). In other words, the tunneling is 
diffusive, if the tunneling probability P(r — r J ,t — t') is 

significant only if |r — r'\/l ~ 1 (where / is a short length 
scale of order of a lattice constant). It should be noted 
that, generally speaking, the question of the conservation 
of in-plane momentum during tunneling is another aspect 
of tunneling not equivalent to the question of coherence 
or incoherence of tunneling (that is the question of the 
dephasing of an electron). Equation || can be rewritten 
in the following form |l(J : 



I(V) 



2eSt 2 , 



J dEdE'dkdk'A^k, E)A 2 (k' ,E')(f{E)(l - f(E'))P(E + eV-E',k- k') - 

-f(E')(l - f(E))P(E' -eV-E,k~ k% 



(3) 



where Ai^(k,E) are the spectral functions, and f(E) is 
a Fermi function. 

The tunneling probability. We need to calculate the 
correlation function: (^e 1 ^^^ e~ 4l ^ r ! * A. The averaging 
can be done if we assume the field ip(r, r) is Gaussian 
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correlated. Using gauge invariance, the phase tp(r, r) can 
be rewritten as ip(r,r) = f^^SV^jtydt, where SV(r, t) 
is the local voltage difference between two planes. This 
way we get an expression for the tunneling probability 
(the calculation is a generalization of the exact calcula- 
tion from Ref.( [1, p.273-277)): 



P{5r 



r'5t 



t-t') 
1, 



exp(-R(6r, St)), 

t-2 



d f 

R{5r, 5t)= j — d 2 q < 8V 2 tU> > (1 - cos{uu8t + $S~r))coth—. 



(4) 



In this paper, we assume that the most effective "de- 
tuning" fluctuations are the voltage fluctuations (or re- 
lated charge fluctuations). It is important to point out 
that the same method can be used to calculate the 
"blocking" of the tunneling due to any mechanism of in- 
plane scattering. Since the nature of the ground state of 
cuprates and therefore the spectrum of the fluctuations is 
not known, later we examine several general forms of the 
spectrum. The problem of incoherent tunneling between 



a couple of two-dimensional planes is a natural general- 
ization of the spin-boson model of quantum dissipation. 

In the view of the importance of the spatial dependence 
of the tunneling probability, we give several different ar- 
guments proving the diffusive nature of the tunneling (if 
i~hop 3> T a b) in the normal state. First of all a qualitative 
argument: if the hopping time is much longer the in-plane 
scattering time, an electron experiences many inelastic 
and elastic scattering processes (both not conserving the 
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direction of the in-plane momentum) before the hopping 
between planes. Thus it is intuitively natural to think 
that the momentum is not conserved after the hopping. 
A straightforward quantitative argument is given by the 
analysis of the function R(8r, St) in the exponent of the 
expression for the tunneling probability (Eqn. ||). To 
separate the spatial and time dependence of the func- 



tion R(Sr, St) , we rewrite the multiplier in the integral 
as (1 — cos(u:St + q x Sr)) = (1 — cos(u>St)) + cos(u>8t)(l — 
cos(q x 5r))—sin(q x 5r)sin(oj5t). The integral with the last 
term vanishes, because this term makes the integral ex- 
pression antisymmetric with respect to integration over 
q x . Thus, we can write 



R(Sr, St) = Ro(Sr, St) + R x (8r, St), 
Ro(St) = J ^ J dq x q y < SV 2 U > coth^(l - cos(w5t)), 

Ri{Sr,St) = 2 f ^ f dq x q y < 5V 2 U > coth^cos(ujSt)sin 2 ((q x 5r) /2). 



(5) 
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The space- independent part Ra(St) is calculated later in 
the paper (see Eqn. ||). Below we calculate the func- 
tion R\ (Sr, St) which describes the spatial dependence of 
the tunneling probability. For this calculation we assume 
that the fluctuations are uncorrelated in two planes (see 
below for a more general discussion), and the interplane 
noise is just twice the in-plane Johnson-Nyquist (JN) 
noise (voltage fluctuations). The Johnson-Nyquist noise 
can be calculated from the spectral density of Coulomb 
noise in the two-dimensional plane valid in the hydrody- 
namic approximation p2|: 



Attctq- 



<T 2 CJ 



4:7T 2 a 2 q 2 



(6) 



where 02 is a two-dimensional conductance, and oq = 
e 2 /Ti. We substitute EqnJ| into the expression (||) for 
R\ (Sr, St) and impose an upper cutoff q c on the q- 
integration to take into account of the fact that the ex- 
pression (Js|) is strictly valid only in the long-wavelength 
limit; thus we take q <C qF (but still of the general order 
of q F )- We also impose on the w-integration a lower cut- 
off u>i, the choice of which will be discussed below. Then, 
taking into account the fact that we are interested in val- 
ues of St which are of the general order of magnitude 
Thop and thus several orders of magnitude larger than 
(c29c) _1 , we see that to a good approximation R\(Sr, St) 
factorizes into a product of a function of Sr and a function 
of St: 



F(Sr) 



Rx{Sr,St) : 
dO / dx- 

JO 



F(5r)G(5t), 
- cos(xcos9) 



™2 J m UJ 



coth (^)cos(ujSt). 



(7) 



The 

q c Sr 



expression F(Sr) is approximately ^(q c Sr) 2 for 
•C 1 and 2irln(q c 5r) for q c 5r ^> 1. As for the 
function G(St), we will see below that the "interesting" 
values of St (i.e. those for which the function Ro(St) 
does not suppress P(Sr, St) too badly) are less or of or- 
der of h~/2irakT, where the dimensionless quantity a is 



typically less or of order of 1. Under these conditions, 
provided htoi <C kT which will be satisfied by our choice 
of ui (cf. below), the integral defining G(5t) is domi- 
nated by its lower limit and approximately given by the 
(^-independent expression 



G(5t) 



2(j Q kT 

7TC72 TlUJl 



(8) 



We will make the choice loi ~ 1/t/ iop , on the grounds that 
once we need to allow for appreciable interplane hopping, 
Eqnjf^ for the noise is no longer applicable (and we ex- 
pect the expression for < SV 2 ^ > to decrease as a higher 
power of lu for ui — > 0, thereby effectively cutting off the 
integral @)). We thus have for Ri(Sr, St) the approxi- 
mate expression 

R^Sr^t) ~ (kTT hop /h){a Q /a 2 )(2/TT)F(q c Sr). 

The salient point, now, is that the quantity kTTh op /n, 
which is essentially the ratio of the ab-plane and c-axis 
conductivities, is of order 10 2 — 10 4 for most of the 
cuprates, while the ratio 02 /o~q is never greater than 
about 10. Thus, the quantity i?i has a value large 
compared to unity for values of Sr small compared to 
l/q c , and we can approximate the expression F(Sr) by 
its limiting form ^(q c 5r) 2 . Thus, the "effective area" 
S e ff = (8r e ff) 2 for which Ri is appreciable is defined by 



Seff 



4(g 2 /o-q) 
q 2 {kTT hop /K) 



and by the above argument this is much smaller than 
1/q 2 and thus at most of the order of l/q F - At distances 
of this order formulas such as (||) should no longer be 
taken seriously, but the crucial upshot of the argument 
is that coherence between tunneling events separated in 
space by more than ~ 1/qF can be simply neglected. To 
put it differently, the tunneling is effectively local (diffu- 
sive) ; the effective rms change in momentum in the course 
of a tunneling event is of order of qp (cf. below). It is 
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noteworthy that this is so even if the momentum cutoff 
q c on the voltage fluctuations is only a small fraction of 
qf ; consideration of shorter- wavelength fluctuations can 
only strengthen this conclusion. 

Another argument estimates directly the change of mo- 
mentum in the process of tunneling. The change of the 
momentum due to a fluctuation of the electromagnetic 
potential is Sp ~ (e8A)/c, therefore < Sp 2 >= (e/c) 2 < 
5 A 2 > . We can relate the correlation function of the vec- 
tor potential with the correlation function of the scalar 
potential by gauge transformation (assuming only lon- 
gitudinal fluctuations): < A 2 ^ >= 2Jr < 5V£ U >. 
Using Eqn. ^ at low frequencies lu < o~2q, the cor- 
relation function can be approximated as < A 2 u >= 

— —— ^— — 5 , g 3 9 ~ 4 ? c . Thus the variance of momentum 

< Sp 2 > is < Sp 2 >=-^jdq qjdui^^ZL. We see 
that the integral over frequency is diverging for oj < T 
as J This integral can be cut off again on 1/t/ iop , 
therefore < Sp 2 >^ <£-^ { T hopkT) (at any finite tem- 
perature). This estimate gives a result equivalent to the 
earlier calculation. This indicates again that the in-plane 
momentum is completely randomized after the process of 
tunneling. All these arguments validate theories of c-axis 
transport in the normal state assuming non-conservation 
of in-plane momentum during tunneling. The fact of the 
non-conservation of momentum fcii in the normal state (if 
i~hop 3> T a b) is quite general, a sufficient condition as can 
be seen from the above discussion is the ohmic density 
of the noise < 5V r (u)) 2 >~ to for ui — > 0. ft is important 
to stress that the "diffusivity" of the tunneling is due to 
specific form of the spectrum of the voltage fluctuations 
(and not due to short links or impurities!), for instance, 
it may not be true in the superconducting state. It is 
important to point out that the question of the in-plane 
momentum conservation for c-axis tunneling can be ex- 
amined experimentally |D| supporting or disproving the 
above arguments. 

When the tunneling is diffusive, the "detuning fluctu- 
ations" are simply the local voltage fluctuations 

< SV 2 >= J d 2 q < SV 2 ^ >= au, 

where a is the microscopic parameter describing the 
ohmic density of the noise. As can be seen below this 
microscopic parameter a is sufficient to describe all dc 
and ac dependencies of the c-axis conductivity. When 
the two planes are widely separated and isolated (the 
situation possibly realized in Bi — 2201), the noise spec- 
tra in each plane are uncorrelated, so that the interplane 
noise spectrum is just the sum of the noise spectra in 
each plane. In such a case, the coefficient a should be 
determined only by the properties of the copper-oxygen 
plane. If the planes are moved closer together, so that 
the inter-plane Coulomb interaction become relevant, the 
intensity of the interplane noise (assuming no interplane 



hopping) increases due to the presence of the acoustic 
(out-of-phase) plasmon in this bi-layer structure. Un- 
fortunately, realistically the noise between planes can be 
suppressed and correlated at low frequencies because of 
interplane hopping and become dependent on the inter- 
layer structure thus implying different values of a for dif- 
ferent cuprate materials. It is known experimentally that 
the detailed temperature dependence of the c-axis resis- 
tivity is very sensitive to several factors (sample prepa- 
ration, interlayer structure and doping). The role of the 
interlayer structure (different intercalating atoms, chains 
and additional layers present in some compounds) and 
the structure of the tunneling matrix element t±(k x , k y ) 
(e.g. in-plane anisotropy) is not clear, it makes the ques- 
tion of the interlayer noise and the c-axis transport prop- 
erties even more complex. 

We have examined several models describing the volt- 
age noise to estimate the microscopic parameter a (see 
the appendix for the discussion of this question, also H|). 
The goal of such exercise is to verify a crude consistency 
between the estimate of a from microscopic noise and the 
parameter a required to describe the c-axis transport de- 
pendencies. The difficulty (not surprising since these cal- 
culations assume Fermi-liquid or diffusive spectra of den- 
sity fluctuations) common to all of the calculations (see 
the appendix) is that the "microscopic" value of a is sig- 
nificantly smaller than the value of order of 1 necessary to 
explain the c-axis transport properties. For our approach 
to be valid a large (a ~ 1) density of "detuning fluctua- 
tions" (not present in RPA or Fermi liquid pictures) is vi- 
tal. An alternative approach is to extract the charge fluc- 
tuation noise directly from the optical reflectivity mea- 
surements. The most dramatic difference between good 
metals and cuprates seen in Raman and optical mea- 
surements of the in-plane dielectric constant e a b(q,u>) is 
that the low-frequency noise for cuprates is ohmic (linear) 
even for q — > . If we write Im(—l/e a b(q,u))) = juj for 
q ~ 0, then from the experimental data pl| 7 ~ 0.2eU _1 . 

Since < SV 2 >= J ^V 2 (q)Im(-l/e k>u ) ~ ^cj, we 
2 

get a ~ %f £ 7- If we take the upper cut-off wave vec- 
tor q c ~ — (a is a in-plane lattice constant), we get 

2 

a ~ —7 ~ 1.2. It shows that a may be of order of 1, 
exactly what is required to explain the c-axis transport 
properties. We now use the "local" approximation justi- 
fied above to calculate the dc and ac c-axis conductivity. 

First of all, we calculate the local tunneling probability 
P(r = 0, t) (or rather its Fourier transform P(e, k — k')). 
The quantity P(e, k—k') can be interpreted as a probabil- 
ity to exchange an energy e and momentum (k — k') with 
fluctuating fields. Since, as argued above, the tunneling 
probability is strongly peaked at r = (if Ty> op 3> T a b), 
the Fourier transform in momentum space P(e, k — k') is 
essentially independent of (k — k'), therefore we omit the 
index (A: — k') below. For the small values of the coef- 
ficient a (a <g; 1) describing the spectral density of the 
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local voltage noise < SV^ >= au>, the tunneling proba- 
bility can be calculated analytically. The function Ro(5t) 
in the exponent of the P(Sr = 0, St) is 

Ro(St) = I ^ < SVl > (1 - cos(uSt))coth(^) = 



30 dtu . co . 

— aw(l - cos{uj5t))coth{-^) 

h 2aln(uj c /kT), 



(9) 



where = l/2wakT and w c is a high-frequency cutoff 
which we associate with the inverse of the transversal 
time U/rtr ~ leV. After a Fourier transform, it gives the 
probability of tunneling P(e) = 5 e// (^) 2a ga+ 2 ( g r)a , 
where Se// is the effective area of tunneling discussed 
above. For a of order of 1, the tunneling probabil- 
ity is strongly suppressed and weakly depends on e (for 
e < kT). It cannot be calculated explicitly analytically, 
but at small e can be approximated as a function inde- 



CJ tl (kT\2a 1 



pendent of e: P(e) 

TTie optical conductivity (j c {lo, kT). The tunneling con- 
ductance o~ c (V). Equation^ can be further transformed 
assuming the diffusive tunneling probability. In this case, 
we can integrate over the in-plane momenta (k,k') sep- 
arately. The next transformation is due to the detailed 



balance condition (see Ref. |11|), eventually we can write 
the expression for the dc conductivity o~ c {V) as a function 
of the applied voltage V: 



= d* ± d,4 D l-e 



-0eV 



V 



do, 

— 00 



From Eqn. [l(J in the limit of small voltage (V <C kT) (for 
all a as long as the tunneling is diffusive) appropriate for 
the dc measurements, we calculate the temperature de- 
pendence of the c-axis conductivity 



<r c (T) = -t ± dv 2D S eff ( — ) — . 



kT\ 



2d 



1 



(11) 



This result describes the c-axis resistivity p c (T) either 
constant or diverging at low temperatures found exper- 
imentally in several compounds (Si-family, Hg — 1201, 
Tl - 2212,17 - 1212 and slightly underdoped La - 214). 

In order to calculate the optical conductivity o~ c (uj), we 
make the following observation. If we apply the external 
dc voltage, the tunneling probability acquires an addi- 
tional phase factor P(t — t') ~ ex p{j- f t i Vdr) (in Eqn||, 
it gives a shift in the energy difference (after a Fourier 
transform) P(E + eV — E')). Thus schematically, the 
conductance is <r c (V) ~ y J dt' e^ v { - t ~ t '} (...). For the ac 
voltage, 

P(t - t') - exvCr f Ve^dr) = 



= exp(-V ; ) W 

n icu 

l + If_(l- e M*-f)). 
n lu 



(12) 



In the linear response (the optical conductivity is a lin- 
ear response) and separating a corresponding harmonic 
of the current proportional to e luJt , we conclude that the 
dependence of the conductivity a c {co) on the frequency 
w is equivalent to the dependence on the voltage V, such 
that a c (u>) = a c (V — > to) (if the tunneling is incoherent 
and diffusive). Namely, 



a c (u>) 



et\dv2 D 1 



-f3u 



+00 



dc 



iP{ia - e) 



1 



(13) 



A ubiquitous nearly flat optical response for cr c (w) is 
observed experimentally (if measurements exist) in these 
compounds. Indeed, a qualitative and numerical analysis 
of the frequency dependence o- c (ui) of Eqn.|l^ (or equiva- 
lently dependence on the voltage) indicates a very weak 
dependence on frequency. 

It is interesting that under the conditions of incoher- 
ent tunneling ohmic I — V curves (constant conductance 
c c (V) as a function of voltage) correspond to a flat op- 
tical conductivity a c (co). The correspondence cf c (uj) = 
a c (V — > lu) can be directly checked experimentally by 
comparing the tunneling conductance <J C (V) = I(V)/V 
and the frequency dependence of o~ c {lo) from the opti- 
cal reflectivity measurements. The graphs of cr c (V) and 
o~ c (lo) can be taken from experimental papers p6[ . It 
appears to be roughly true for optimally doped and un- 
derdoped compounds in the normal state. 

Discussion of experiments. The above results are in- 
sufficient to describe all experimental data for different 
dopings. In optimally doped La — 214 and YBaCuO and 
other overdoped cuprates the c-axis resistivity has a lin- 
ear temperature dependence with a large "residual" value 
(intercept at T = OK) [0. In these compounds the hop- 
ping time becomes comparable to the in-plane scattering 
time. Thus we can assume that specular tunneling be- 
comes possible. In this crossover situation between two 
limiting pictures of the tunneling between planes and the 
anisotropic band along c-axis, we can think about two 
channels of conduction. One channel is diffusive, while 
another one is specular. The tunneling probability is the 
sum of probabilities to tunnel without and with the con- 
servation of the in-plane momentum. In this case, the 
total c-axis conductivity is the sum of conductivities in 
each channel. If the fermi surface has very anisotropic 
properties, electrons from one part of the fermi surface 
can tunnel specularly, while electrons from other parts 
tunnel diffusively. In this scenario, the second channel 
of conduction (conserving fen) can be due to the diag- 
onal parts of the fermi-surface in the normal state. It 
appears from photoemission experiments that the quasi- 
particles along diagonals of the Brillouin zone {k x — zbfc^) 



G 



have longer life-times T a b,diag (which should be compared 
with the hopping time). These quasiparticles can tunnel 
then with conserved momentum. It is natural to suggest 
that for overdoped cuprates the c-axis transport is dom- 
inated by incoherent, but specular channel, while for un- 
derdoped cuprates the diffusive channel is only present. 

We can calculate the conductivity of a specular chan- 
nel, if we assume that the time dependence of the spec- 
ular tunneling probability is P{5t) ~ exp{~^) with 
calculated for a weak detuning, that is = l/2nakT. 
Thus we substitute the tunneling probability of the form 
P 2 (e - e', k - k') = S(k - k') {e ^,f^ akT) , (or any form 

P 2 (e-e',k-k') ~ 5{k- k')^f(^f)) to the Eqn. § we 
get the contribution to the c-axis conductivity 



e 2 A 
<T2,c( T ) = T^^Pf' 



(14) 



where A is a numerical coefficient. This result is a well- 
known result for incoherent tunneling with conservation 
of the in-plane momentum |jq| . It is suggested to ex- 
plain the linear temperature dependence of c-axis resis- 
tivity observed in some compounds P9]. Another im- 
portant consideration is that the band calculations pre- 
dict a significant angular dependence of t± in some fam- 
ilies of cuprates ||. Due to this reason, the contribution 
from specular channel (from diagonal parts of the fermi- 
surface) can be reduced. 

Yet another complexity of cuprates with multiple 
planes per unit cell is the question of intra-unit cell and 
inter-cell conduction. In this case, the resistance as- 
sociated with hopping between planes of the unit cell 
(intra-cell resistance) and the resistance associated with 
hopping between different unit cells (inter-cell resistance) 
should be discussed. The total resistance is certainly a 
sum of intra- and inter- unit cell resistances. It may not be 
correct to assume (as is frequently done in the literature) 
that the intra-cell resistance is negligible; a systematic 
experimental investigation of this question is necessary. 
We hope to discuss this question elsewhere. 

In conclusion, a picture of c-axis interlayer (and inter- 
cell) tunneling strongly suppressed by voltage fluctua- 
tions is proposed. This approach can provide a consistent 
understanding of observed temperature and frequency 
dependencies of c-axis conductivity in the normal state. 
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Appendix. 

In this appendix, we demonstrate an example of a cal- 
culation of the parameter a describing the low frequency 
voltage fluctuations based on a particular microscopic 
model of the density fluctuation spectrum. At the end of 
the appendix, we list several results calculated for various 



other microscopic models. 

If the two planes are widely separated, then the voltage 
fluctuations can be treated independently in each plane. 
In this case, the interplane noise is just twice the in- 
plane local Johnson-Nyquist noise in each plane. For a 
high-density electron gas in a hydrodynamic approxima- 
tion (cjt <C 1 and ql <C 1), the charge density-density 
susceptibility of the two-dimensional electron gas can be 
written as |20|] 



X(k,u>) 



s 2 kTFk 2 



1 



27re 2 lo(lu + i/r) — s 2 k 2 — s 2 kTFk ' 



where s 2 = v 2 F /2, kf p — 2 ™°i is the Thomas- Fermi 
wave number, and t is a phcnomcnological relaxation 
time. The spectral density of the voltage fluctuations is 



<VL >=V 2 Imx(k,Lu), 



where Vi- 



is the two-dimensional Coulomb in- 



27re- 

teraction. Eventually, we need to calculate the partial 
frequency-dependent spectral density of the voltage noise 
< V 2 >= J ^ < Vg u >. The calculation gives 



< Kf >= 



dkk (2Tre 2 ) 2 s 2 k TF k 
IhT fc 2 27re 2 
W (l/r) 



{lo/t) 2 + (s 2 k 2 + s 2 k T Fk) 2 

a 2 „2^,2 



), 



- 2 -u)ln{- 

where uj c is an infra-red cut-off frequency. It implies to 
the accuracy of the value of the logarithm that the pa- 
rameter a is 



sUtftU 2€ft 
We can rewrite this expression in the following form: 

a = 

7T 2 D 



is the 



where a Q = e 2 /h, o 2 d = e 2 v 2 DD = | e ( ^ t) 
two-dimensional conductivity (it can be shown to have 
a Drude frequency dependence). The above results may 
be used to estimate the value of the parameter a in the 
cuprates. Two difficulties can be seen from such literal 
application of the above model. First of all, the value 
of a is significantly smaller than one. Second, and more 
importantly, the parameter a appears to be temperature 
dependent. It should be realized that the spectrum of 
charge fluctuations in the cuprates is much more com- 
plex and not represented correctly by the above simple 
model. We investigated several other simple microscopic 
models in order to get further insight into this question. 
At this moment, it seems more reasonable to extract the 
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charge fluctuation spectrum directly from experiment as 
shown in the text of this paper. 

The results for several other microscopic models are 
summarized below. If the spectrum of charge fluctuations 
is given by the weakly damped acoustic two-dimensional 



plasmon, then the parameter a is a 



where 



s 2 (l + k TF d) ' 

d is the inter-plane distance. Another calculation tak- 
ing the voltage noise due to electron-hole pairs of the 
two-dimensional Fermi liquids gives a ~ 1/Att. If the 
interplane Coulomb interaction is taken into account in 
the RPA approximation for the same calculation (el-hole 
pairs), a — l/(27rfcjrd). We hope to present the details 
of calculations and expanded arguments elsewhere. 
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